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^ ; A new thermalization scenario for heavy ion collisions is discussed. It is based on 
the Hawking-Unruh effect: an observer moving with an acceleration a experiences the 
influence of a thermal bath with an effective temperature T = a/27r, similar to the one 
present in the vicinity of a black hole horizon. In the case of heavy ion collisions, the 
acceleration is caused by a pulse of chromo-electric field E ~ Q 2 S / g (Q s is the saturation 
scale, and g is the strong coupling), the typical acceleration is a ~ Q s , and the heat bath 
temperature is T ~ Q s /2tt ~ 200 MeV. In nuclear collisions at sufficiently high energies 
the effect can induce a rapid thermalization over the time period of r ~ tt/Q s accompanied 
by phase transitions. A specific example of chiral symmetry restoration induced by the 
chromo-electric field is considered; it is mathematically analogous to the phase transition 
occurring in the vicinity of a black hole. 



1. INTRODUCTION 

There is something puzzling in the data on hadron production in various processes, 
from e + e~ annihilation and deep-inelastic scattering to heavy ion collisions: the relative 
abundances of different hadron species appear to follow the statistical distribution with a 
surprising accuracy (for reviews, see e.g. [HJI2])- Moreover, at small transverse momenta 
the spectra of the produced hadrons also look approximately thermal. While in heavy 
ion collisions it is possible to expect the emergence of statistical distributions as a result 
of intense re-interactions between the produced particles, this seems very implausible in 
e + e~ annihilation at high energies where the process of hadronization is stretched 
in space over a long distance ~ and the density of produced hadrons is small 

(/i ~ Aqcd is an infrared cut-off describing the hadronization scale). 

Statistical distributions may emerge as a result of the saddle-point approximation to 
the multi-particle phase space, when the dynamics is inessential and the production 
mechanism is "phase space dominated" (see, for example, [[3] and references therein). 
However, in e + e~ annihilation this mechanism can hardly be expected to work: the jet 
structure, the angular distributions of the produced hadrons, and inter-jet correlations 
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point to the all-important role of QCD dynamics of gluon radiation (for a recent review, 
see [E]), and thus the "phase space dominance" cannot be invoked. It thus looks natural 
to think that the emergence of statistical hadron abundances has something to do with 
the process of hadronization - in other words, with the way in which the QCD vacuum 
responds to the external color fields, as was advocated by Dokshitzer [ |1] . 

This attractive and deep idea however does not yet allow one to understand why the 
vacuum response to the external carriers of color field leads to the emergence of statistical 
distributions. To do this, we are forced to discuss the structure of the QCD vacuum. Little 
is known about it, but we do know that QCD vacuum is populated by the fluctuations 
of color fields, some of which are of semi-classical nature [E]. It may therefore be useful 
to examine what is known at present about the interactions of charged quanta with 
background classical fields. 

This logic will first lead us to the discussion of quantum fluctuations in the background 
of the gravitational field of a black hole, where the quantum radiation appears to be 
thermal as shown by Hawking [EJ- It was demonstrated by Unruh [|7j that the Hawking 
phenomenon should be present in any non-inertial frame; indeed, this is required by 
Einstein's Equivalence Principle. Then we will proceed to the process of electron-positron 
pair production in the background electric field, analyzed in QED by Schwinger [ 8 . We 
will see that the Schwinger formula for the rate of e + e~ pair production in a constant 
electric field allows for a simple statistical interpretation; moreover, for the case of a 
time-dependent field pulse, the spectra of the produced particles become thermal. 

These two examples inspire the picture in which the thermal character of hadron pro- 
duction emerges from the interactions with the vacuum chromo-electric field, the strength 
of which in the color flux model is parameterized in terms of the string tension, or the 
slope of Regge trajectories b. Indeed, we will see how the well-known formula for the 
Hagedorn temperature 

J- Hagedorn ^ ^ 

can be derived in this way. 

What happens if we create semi-classical color fields of strength exceeding the strength 
of the vacuum fields? This can be achieved in the collisions of heavy ions at high energies, 
which are accompanied by a short pulse of the chromo-electric field E ~ Q 2 s /g of duration 
r ~ Qj 1 ] here Q s is the saturation scale Q s in the Color Glass Condensate (for a review at 
this Conference, see [El), an d g is the strong coupling. I will argue that such a strong color 
field induces the creation of partons with a distribution which is isotropic and thermal in 
a co-moving local frame, with an effective temperature T in ~ Q s /2tt. For high enough 
energies and heavy enough nuclei, the value of T in exceeds the Hagedorn temperature; the 
produced thermal matter is thus in the deconfined phase. The phase transition in this 
case can also be understood in a geometrical picture, in which the acceleration a ~ Q s 
in the chromo-electric field determines the curvature of space in the non-inertial Rindler 
frame. This is mathematically analogous to the phase transition induced by the presence 
of a massive black hole. This talk is based on the paper with Kirill Tuchin [Ell, to which 
I refer for details and references. 
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2. QUANTUM FLUCTUATIONS IN THE CLASSICAL BACKGROUND 

2.1. Black hole evaporation 

In 1974 Hawking [ ,5J demonstrated that black holes evaporate by quantum pair pro- 
duction, and behave as if they have an effective temperature of 

T„ - £, (2) 

where k = 4M is the acceleration of gravity at the surface of a black hole of mass M. 
The thermal character of the black hole radiation stems from the presence of the event 
horizon, which hides the interior of the black hole from an outside observer. The rate 
of pair production in the gravitational background of a black hole can be evaluated by 
considering the tunneling through the event horizon. Parikh and Wilczek [ showed 
that the imaginary part of the action for this classically forbidden process corresponds to 
the exponent of the Boltzmann factor describing the thermal emission 2 . 

Unruh [Cj has found that a similar effect arises in a uniformly accelerated frame, where 
an observer detects an apparent thermal radiation with the temperature 

Tu = h (3) 

(a is the acceleration). The event horizon in this case emerges due to the existence of 
causally disconnected regions of space-time, conveniently described by using the Rindler 
coordinates. 



2.2. Pair production in a constant electric field 

The effects associated with a heat bath of temperature (J3J) usually are not easy to 
detect because of the smallness of the acceleration a in realistic experimental conditions. 
For example, for the acceleration of gravity on the surface of Earth g ~ 9.8 m s~ 2 the 
corresponding temperature is only T ~ 4 x 10~ 20 K. (The wavelength of the thermal 
radiation in this case is about 1 parsec - so it would have to be detected on a cosmic 
scale!) 

Much larger accelerations can be achieved in electromagnetic fields, and Bell and 
Leinaas [^j considered the possible manifestations of the Hawking-Unruh effect in parti- 
cle accelerators. They argued that the presence of an apparent heat bath can cause beam 
depolarization. Indeed, if the energies of spin-up and spin-down E± states of a particle 
in the magnetic field of an accelerator differ by AE = E^ — E^, the Hawking-Unruh effect 
would lead to the thermal ratio of the occupation probabilities 

_.ex P (-— ), (4) 

where the Unruh temperature is determined by the particle acceleration. According 
to (J3J), a pure polarization state of a particle in the accelerator is inevitably diluted by 
the acceleration. 



2 Conservation laws also imply a non-thermal correction to the emission rate [ lllj. possibly causing a 
leakage of information from the black hole. 
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If the energy spectrum of an accelerated observer is continuous, as is the case for 
a particle of mass m with a transverse (with respect to the direction of acceleration) 
momentum p?, a straightforward extension of (jlj) leads to a thermal distribution in the 



"transverse mass" rriT = ym 2 + Pt- 

W m (p T ) ~ exp (-jPj ■ (5) 

An important example is provided by the dynamics of charged particles in external electric 
fields. Consider a particle of mass m and charge e in an external electric field of strength 
E. Under the influence of the Lorentz force, it moves with an acceleration a = eE/m, and 
the corresponding temperature is Tjj = a/2n. The Boltzman factor exp(— m/Ty) entering 
the particle creation rate in this case is 

W^-expj- — —J. (6) 

This expression looks familiar - in fact, © differs from the classic Schwinger result for 
the rate of particle production in a constant electric field only by a factor of two in the 
exponent. 

Is this a coincidence? To answer this question, let us have a closer look at the corre- 
sponding action: 

S — J ( — mds — ecpdt) , (7) 

where (p is the electric potential, which for a constant electric field aligned along the x 
axis is <p = —Ex modulo an additive constant; the invariant interval is given by ds 2 = 
(1 — v 2 {t)) dt 2 . Using the equations of motion, we can evaluate the action to find 

S(t) = ^ dt(-myjl - v(t) 2 + eExit)) 

= arcsinh(ar) + — - (ar(Vl + a 2 r 2 — 2) + arcsinh(ar)) + const. (8) 

a 2 a 2 v 7 

In classical mechanics the equations of motion completely specify the trajectory of a 
uniformly accelerating particle moving under the influence of a constant force F = —eWtp. 
In contrast, in quantum theory the particle has a finite probability to be found under the 
potential barrier V(x) = eEx in the classically forbidden region. Mathematically, it comes 
about since the action (JHJ) being an analytic function of r has an imaginary part in the 
Euclidean space 

_ . mix eEn urn 2 

ImSM = — -^=271' < 9) 

The imaginary part of the action (jSJ) corresponds to the motion of a particle in Euclidean 
space along the trajectory 



x(t E ) = a" 1 (a/1 - a 2 t\ - 1). (10) 
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Note that unlike in Minkowski space the Euclidean trajectory is bouncing between the 
two identical points x a = —a^ 1 at t S a = — a -1 and Xb = — a -1 at t^,b — an d the 
turning point x a = at t E)CL = 0. Using (jSJ) we can find the Euclidean action between the 
points a and 6; it is given by SE(x(t E )) = ImS = nm 2 /2eE. 

It is well known that a quasi-classical exponent describing the decay of a metastable 
state is given by the Euclidean action of the bouncing solution, (fTTj) . The rate of tunneling 
under the potential barrier in the quasi-classical approximation is thus given by 



W£ ~ exp(-2ImS) = exp j . (11) 

Equation (JTTJl gives the probability to produce a particle and its antiparticle (each of mass 
m) out of the vacuum by a constant electric field E; note that the incorrect factor of 2 in 
the exponent of (JHJ) has now disappeared due to the contribution of the field term in the 
action (jJJ. 

The ratio of the probabilities to produce states of masses M and m is then 

The relation (|T2*Jl allows a dual interpretation in terms of both Unruh and Schwinger 
effects (see e.g. [ H21 EH EB UHj and references therein). First, consider a detector with 
quantum levels m and M moving in a constant electric field. Each level is accelerated 
differently, however if the splitting is not large, M-ra<mwe can introduce the average 
acceleration of the detector 

- (13) 



M + m 

Substituting (JT3j) into (|T2*|) we arrive at 
W§ [2tt(M - m) 



exp ^- '- . (14) 



This expression is reminiscent of the Boltzmann weight in a heat bath with an effective 
temperature 0: T = a/2ir. It implies that the detector is effectively immersed in a 
photon heat bath at temperature T w eEjuvn. This is the renown Unruh effect [|7]. 

It is important to remember that for a constant electric field the momentum distribution 
of the produced charged particles allows a statistical interpretation only in the transverse 
to the field direction. However for a short pulse of an electric field of duration r < m/ (eE) 
the distribution becomes thermal in all three directions (see [EJ and references therein) - 
physically, this happens because the field has not enough time to perform work on curving 
the momenta of the produced particles. 

Let us now discuss ()12 |) .([TT j) from the viewpoint of a field-theoretical derivation done 
by Schwinger [ |H] . It is clear that ([lip cannot be expanded neither in powers of the 
coupling e, nor in powers of the field E, and so cannot be reproduced in any finite order 
of perturbation theory. Schwinger considered the one-loop QED action describing the 
electron-positron fluctuations in the background of the external electric field E. He has 
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found that the series in the number of external field insertions indeed diverges (is not 
Borel summable). As a result the action ceases to be real, and develops an imaginary 
part, similarly to the simple semi-classical example considered above. 

The Hawking-Unruh interpretation therefore appears to capture an essentially non- 
perturbative dynamics. Indeed, a more rigorous treatment (for a review, see e.g. [ I16j) 
allows to establish that the Bogoliubov transformations relating the particle creation and 
annihilation operators in Minkowski and Rindler spaces describe a rearrangement of the 
vacuum structure which cannot be captured by perturbative series. 

3. Event horizon and thermalization in high energy hadronic interactions 

3.1. Unruh effect, Hagedorn temperature, and parton saturation 

We are now ready to address the case of hadronic interactions at high energies, which is 
the main subject of this talk. Consider a high-energy hadron of mass m and momentum 
P which interacts with an external field (e.g., another hadron) and transforms into a 
hadronic final state of invariant mass M 3> m. This transformation is accompanied by a 
change in the longitudinal momentum 

q L = \JE 2 - m 2 - VE 2 - M 2 ~ — , (15) 

and therefore by a deceleration; we assumed that the particle m is relativistic, with energy 
E ~ p. 

The probability for a transition to a state with an invariant mass M is given by 

P(M <- m) = 2tt|T(M <- m)| 2 p(M), (16) 

where T(M <— m) is the transition amplitude, and p(M) is the density of hadronic final 
states. According to the results of the previous section, we expect that under the influence 
of deceleration a which accompanies the change of momentum (|15jl. the probability |T| 2 
will be determined by the Unruh effect and will be given by 

|T(M <- m)| 2 ~ exp(-27rM/a) (17) 

in the absence of any dynamical correlations; we assume M ^> m. 

To evaluate the density of states p(M), let us first use the dual resonance model (see 
e.g. [UH, [EE]), in which 

p(M) ~ exp ( i| Vb Mj , (18) 

where b is the universal slope of the Regge trajectories, related to the string tension a by 
the relation a = 1/ (27r6). 

The unitarity dictates that the sum of the probabilities ([16)1 over all finite states M 
should be finite. Therefore, by converting the sum into integral over M one can see that 
the eqs ()17j) and (|T5|) impose the following bound on the value of acceleration a: 
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The quantity on the r.h.s. of (|19|) is known as the Hagedorn temperature [EI! ~~ the "limit- 
ing temperature of hadronic matter" derived traditionally from hadron thermodynamics. 
In our case it stems from the existence of a "limiting acceleration" qq: 




°>o = \h; b f - (20) 



The meaning of the "limiting temperature" in hadron thermodynamics is well-known: 
above it, hadronic matter undergoes a phase transition into the deconfmed phase, in 
which the quarks and gluons become the dynamical degrees of freedom. To establish the 
meaning of the limiting acceleration (|20J) . let us consider a dissociation of the incident 
hadron into a large number n ^> 1 of partons. In this case the phase space density (jl8)l can 
be evaluated by the saddle point method ("statistical approximation"), with the result 
(see e.g. [ED|) 

p(M) ~ exp(/3M); (21) 

where (3~ l is determined by a typical parton momentum in the center-of-mass frame of 
the partonic configuration. When interpreted in partonic language, eq (jl8|) thus implies a 
constant value of mean transverse momentum px ~ ~ b^ 1 ^ 2 . On the other hand, in 
the parton saturation picture, the mean transverse momentum has to be associated with 
the "saturation scale" Q s determined by the density of partons in the transverse plane 
within the wave function of the incident hadron (or a nucleus). This leads to the phase 
space density p(M) ~ exp(M/Q s ). The unitarity condition and the formulae (|T7| . (fTKj) 
thus lead us to the acceleration a = Q s , which can exceed ()2U|) . and to the conclusion that 
the final partonic states are described by a thermal distribution with the temperature 

T = — . (22) 

The same result can be obtained by considering the acceleration a = gE/m of a parton 
with off-shellness m = y/p 1 ~ Q s in an external color field gE — Q 2 S . It is interesting 
to note that to exceed the limiting acceleration ([20)1 . and thus the limiting Hagedorn 
temperature (j!9|) for the produced hadronic matter, one has to build up strong color 
fields, exceeding gE ~ 1/6. This is achieved by parton saturation in the Color Glass 
Condensate, when gE ~ Q 2 S > gE Q at sufficiently high energies and/or large mass numbers 
of the colliding nuclei. Parton saturation in the initial wave functions thus seems to be a 
necessary pre-requisite for the emergence of a thermal deconfined partonic matter in the 
final state. 

The thermal distribution is built over the time period of 

^therm ~ T 1 = — . (23) 
^% s 

As discussed above, this apparent thermalization originates from the presence of the event 
horizon in an accelerating frame: the incident hadron decelerates in an external color field, 
which causes the emergence of the causal horizon. Quantum tunneling through this event 
horizon then produces a thermal final state of partons, in complete analogy with the 
thermal character of quantum radiation from black holes. 
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3.2. The space-time picture of relativistic heavy ion collisions 

The conventional space-time picture of a relativistic heavy ion collision is depicted in 
the left panel of Fig^ The colliding heavy ions approach the interaction region along 
the light cone from x = t = — oo and x = — t = oo. The partons inside the nuclei in the 
spirit of the collinear factorization approach are assumed to have a vanishing transverse 
momentum fcy, have a zero virtuality k 2 = —k\ = 0, and thus are also localized on the 
light cone at ±x = t. The collision at x — t — produces the final state particles with 
transverse momenta px which according to the uncertainty principle approach their mass 
shell at a proper time r = (t 2 — x 2 ) 1 ^ 2 = r ~ 1/pt- 




Figure 1. The space-time picture of a relativistic heavy ion collision. Left: heavy ions 
approach the interaction region around x = t = along the light cone from x = t = — oo 
and x = — t = oo. The collision at x = t — produces the final state particles which 
approach the mass shell at some proper time r = (t 2 — x 2 ) 1 ^ 2 (or, equivalently, along the 
surface of Rindler space p 2 = pi = x 2 — t 2 < 0). The produced particles are distributed 
in rapidity, or in the Rindler coordinate r\ — |ln |^ . Right: also shown in the left (L) 
and right (R) sectors are the trajectories of space-like p 2 ~ — Q 2 < partons confined 
in the initial nuclear wave functions characterized by the saturation scale Q s . Quantum 
tunneling of partons from left (L) and right (R) sectors through this event horizon into 
the future (F) indicated by dashed arrows creates a thermal state of parton matter with 
the temperature T = Q s /2ti. 



For further discussion it is convenient to introduce the Rindler coordinates 



P 



x 



t 1 



V 



In 



t + x 



t 



X 



(24) 



The surface of a fixed proper time which is a hyperbola in Minkowski space thus represents 
a line at p 2 = x 2 — t 2 = pjj < in Rindler space. The Rindler coordinate r] in high energy 
physics is often called a space-time rapidity. 
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Consider now the case when partons in the wave functions of the colliding nuclei have 
non- vanishing transverse momenta, as in the Color Glass Condensate picture where their 
transverse momenta are on the order of the saturation scale Q s . In this case the partons 
are space-like k 2 = —k\ and are located off the light cone. As the colliding nuclei 
approach each other, these partons begin to interact; note that since they are space-like, 
their interactions are acausal, and are responsible for the breakdown of factorization for 
the parton modes with kj> < Q s - The interactions of partons with the color fields gE ~ Q 2 
of another nucleus decelerate them, with a typical acceleration |a| ~ Q s . The space-time 
trajectories of the partons are thus given by hyperbolae in Minkowski space, or by the 
lines with a fixed value of p 2 = — p§ = Q 2 > in Rindler space. These trajectories 
are shown on the right panel of FigO For partons moving in the left (L) and right 
(R) sectors of space-time with an acceleration \a\ = p~ l the light cone surfaces p 2 = 0, 
7] = ±oo represent the event horizon of the future (F). The information from the future 
is hidden from them, and the sector F is classically disconnected from L and R. However, 
as discussed above, the future sector F can be reached from the left L and right R sectors 
by quantum tunneling. 

3.3. Hawking— Unruh phenomenon in the parton language 

It is important to mention that FigC]does not imply any "bouncing" of the colliding 
nuclei - the on-shell particles, as well as the high-momentum valence partons from the 
colliding nuclei are transparent for each other. The soft components of the parton wave 
functions however do interact strongly. Nevertheless, the whole picture at first glance 
looks completely orthogonal to the conventional parton model. Indeed, in parton model 
the Weizsacker- Williams gluon fields surrounding the valence quarks are transverse, with 
E _L H _L p (p is the momentum of the quark). The corresponding "equivalent gluons" 
are almost on mass shell, and no longitudinal chromo-electric field is present: the gluon 
field tensor is flat in the longitudinal direction: F + _ = ("+" and "— " refer to the 
light-cone components). 

However, a more careful analysis reveals that this picture is not complete: the configu- 
ration of the gluon field produced in the sector "F" of FigHJis characterized by F + _ ^ 0, 
with a substantial longitudinal chromo-electric field E z . In the Color Glass Condensate 
picture, the strength of the field is E z ~ Q 2 /g, and the duration of the pulse is ~ Qj 1 . 
Note that in the conventional string model picture, the produced chromo-electric field is 
purely longitudinal, with the strength proportional to the string tension. This exhibits a 
possible continuity between the string and parton approaches to multi-particle produc- 
tion, and suggests the existence of the minimal allowed value of the saturation momentum 
Qs min- Basing on the arguments given above and (JIHJ), one is led to the conclusion that 

Qs min ^^T-Hagedarn — 1 GeV. (25) 

A quantitative analysis of parton production in the longitudinal chromo-electric field 
performed in the framework of the Color Glass Condensate approach is underway [I21j. 

4. SUMMARY 

In this talk I have argued that the statistical features of multi-particle production may 
emerge as a consequence of the Hawking-Unruh effect. The acceleration, and the emer- 
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gence of the corresponding event horizon for partons, is caused by the pulse of chromo- 
electric field which accompanies inelastic interactions at high energies. For hadron inter- 
actions at moderate energies, the effective temperature appears equal to the Hagedorn 
temperature (|I5|). Once the strength of the chromo-electric field E ~ Q 2 s /g exceeds a 
critical value determined by (|25|). the partons are produced with an effective temperature 
T > T Hagedorn , i.e. in a deconfined state. The subsequent evolution of the produced 
partonic system has to be taken into account in this case. 

I am grateful to K. Tuchin for sharing the fun of thinking about the problem discussed 
here, and to T. Csorgo, G. Dunne, R. Glauber, E. Levin, L. McLerran, G. Nayak and R. 
Venugopalan for helpful discussions. 
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